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Abstract. The NA60 experiment has measured the production of muon pairs 
and of charged particles in In+In collisions at a beam energy of 158 A GeV. For 
invariant dimuon masses below the cj> the space-time averaged p spectral function 
was isolated by a novel procedure. It shows a strong broadening but essentially no 
shift in mass. The production of 3/ip was measured as a function of the collision 
centrality. As in previous experiments studying Pb+Pb collisions an anomalous 
supression is observed, setting in at approximately 90 participant nucleons. Using 
the charged particles the reaction plane was reconstructed. The elliptic flow of 
charged particles increases with pt showing a saturation for pt > 2 GeV/c. For 
the first time azimuthal distributions for J/ip are shown. 



PACS numbers: 25.75.-q,25.75.Dw,25.75.Ld,25.75.Nq 

1. Introduction 

High-energy heavy-ion collisions provide a unique opportunity to study the behaviour 
of strongly interacting matter at high densities and at high energy densities. They 
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are the only way to investigate a possible phase transition from hadronic matter to 
a plasma of deconfmed quarks and gluons as well as the restoration of the chiral 
symmetry which is spontaneously broken in the hadronic world. These issues have 
been studied by several experiments at the SPS (Super Proton Synchrotron) at CERN 
and many exciting results were obtained. NA60 is a second generation experiment 
designed to answer specific questions that still remain open after the completion of 
these previous experiments. 

The main focus of NA60 is the study of muon pair production. The experimental 
setup is based on the muon spectrometer previously used by NA38 and NA50 P, 
separated from the target region by a 5.5 m long hadron absorber (mostly carbon) and 
composed of a toroidal magnet, of eight multi-wire proportional chambers and of four 
scintillator trigger telescopes. A Zero-Degree Calorimeter (ZDC) measures the energy 
of the spectator nucleons to determine the collision centrality. The muon spectrometer 
is complemented by a high-granularity silicon pixel tracker of unprecedented radiation 
tolerance [21 13 which was constructed using ALICE sensors and readout chips @|. The 
tracker, embedded in a 2.5 T dipolc magnet in the vertex region, tracks all charged 
particles before the hadron absorber and determines their momenta independently of 
the muon spectrometer. It allows for a determination of the primary interaction vertex 
with an accuracy of 10 — 15 /im in the transverse plane and 200 /xm along the beam 
axis. The matching of the muon tracks before and after the hadron absorber, both in 
angular and in momentum space, improves the dimuon mass resolution in the region 
of the vector mesons uj and <j> to ~20 MeV/c 2 , significantly reduces the combinatorial 
background due to ir and K decays and makes it possible to measure the muon offset 
with respect to the interaction vertex 0. Measuring all charged particles within the 
vertex tracker acceptance allows in addition for the reconstruction of the direction of 
the impact parameter vector and hence to determine the orientation of the reaction 
plane and to study azimuthal angle distributions of the emitted particles. 

The experiment has taken data in 2003 (In+In collisions) and in 2004 
(proton+nucleus collisions). The results reported upon in this paper were obtained 
from the analysis of the data taken in 2003 with a 158 AGeV Indium beam, incident 
on a segmented Indium target of seven disks with a total of 18% nuclear interaction 
length. At an average beam intensity of 5-10 7 ions per 5 s burst, about 3 • 10 12 ions 
were delivered to the experiment and a total of 230 • 10 6 dimuon triggers were recorded. 

2. The p spectral function 

In the low-mass sector (m < 7ti^) the CERES experiment has studied the production 
of electron pairs in p+Be/Au, S+Au and Pb+Au collisions [H]. In nuclear collisions a 
clear excess above the expected hadronic sources has been observed in the dielectron 
mass distribution. The origin of the excess has been commonly interpreted in terms 
of thermal production from the dense hadronic gas created in the collision, mainly 
occurring via the it + tt~ —tp—t c + c~ process. Even if it is clear that an in-medium 
modification of the p must be introduced in order to explain the results, the lack 
of statistics and mass resolution have prevented any detailed understanding of the 
character of the in-medium changes. 

Figure^shows the opposite-sign dimuon mass distribution as measured by NA60, 
integrated over all collision ccntralities. The combinatorial background of uncorrelatcd 
muon pairs mainly originating from 7r and K decays is determined using a mixed-event 
technique 0- After subtraction the remaining opposite-sign muon pairs still contain 
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"signal" fake matches (associations of muons measured in the muon spectrometer 
to non-muon tracks in the vertex tracker), a contribution which is only 7% of the 
combinatorial background level. It has been determined by an overlay Monte Carlo 
method. After subtracting the combinatorial background and the signal fake matches, 
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Figure 1. Mass distribution of opposite-sign dimuons (upper histogram), 
combinatorial background (dashed), signal fake matches (dashed-dotted), and 
resulting signal (lower histogram with error bars). 

the resulting net spectrum contains about 360 000 muon pairs in the mass range 
of the figure, roughly 50% of the total available statistics. The analysis is done in 
four classes of collision centrality defined by the charged particle multiplicity density: 
peripheral (4-30), semipcriphcral (30-110), semiccntral (110-170) and central (170- 
240). Numerically, the average rapidity density in each class is, within 10%, equal 
to the average number of participants in the class. The signal-to-background ratios 
associated with these classes are 2, 1/3, 1/8 and 1/11, respectively. 

The most peripheral data can be described on the basis of known sources, i.e. the 
two-muon and the Dalitz decays of the various vector mesons, plus a contribution from 
open charm decays. The expected mass shape for the various physics processes has 
been obtained propagating decay muons through the NA60 set-up, using GENESIS [7] 
(originally based on [5]) as the event generator and GEANT for tracking. This 
'hadronic' cocktail has been directly fitted to the data and is found to reproduce 
the observed spectrum quite well (Hj- 

In the more central bins, the data can no longer be described on the basis of the 
standard hadron decay cocktail alone, but are indicative of the existence of an excess 
yield. Since the particle ratios are expected to be different from the peripheral data, 
global fits to the more central data are bound to bias both the extracted cocktail 
parameters and an excess with a priori unknown characteristics. We have therefore 
used a novel procedure which is illustrated in Fig. Eland described in detail in |10j . 
The excess is isolated by subtracting the cocktail, without the p, from the data. The 
cocktail is fixed, separately for the major sources and in each centrality bin, by a 
"conservative" approach. The yields of the narrow vector mesons u and <f> are fixed 
so as to get, after subtraction, a smooth underlying continuum. For the 77, an upper 
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Figure 2. Isolation of an 
excess above the electromag- 
netic decays of a "conserva- 
tive" hadron decay cocktail 
(see text). Total data (open 
circles), individual cocktail 
sources (solid), difference 
data (thick triangles), sum 
of cocktail sources and differ- 
ence data (dashed). 
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Figure 3. Comparison of 
the excess mass spectrum for 
the semi-central bin to model 
predictions, made for In-In 
at d7V ch /dr/=f40. Cocktail 
p (thin solid), unmodified p 
(dashed), in-medium broad- 
ening p [n][I2] (thick solid), 
in-medium moving p related 
to Q21 EH (dashed-dotted). 



limit is defined by "saturating" the measured data in the region close to 0.2 GeV/c 2 ; 
this implies the excess to vanish at very low mass by construction. 

The excess mass spectrum resulting from the subtraction of the "conservative" 
hadron decay cocktail from the measured data is shown in Fig. [3] for semi-central 
collisions. The qualitative feature of the spectrum is striking: a peaked structure is 
seen around the position of the nominal p pole. The excess spectrum is consistent 
with an interpretation of the excess as being dominated by rnr annihilation. 

Figure as well shows predictions from theoretical calculations employing the 
broadening scenario of ^2 E] an d the moving-mass scenario related to ^3 I14j . 
Both are evaluated for In+In at dN^/dij = 140 within the same fireball evolution. 
In addition the unmodified p is shown. The data shown in this paper have not 
been corrected for the mass- and pt-dependent acceptance of the NA60 setup. The 
theoretical calculations were therefore propagated through the acceptance filter to 
allow for a fair comparison with the data. The integrals of the theoretical spectra are 
normalized to the data in the mass interval M < 0.9 GeV/c 2 . The unmodified p is 
clearly ruled out. The specific moving-mass scenario plotted here, which fitted the 
CERES data [HI El Ej ; is ruled out as well. The broadening scenario appears more 
realistic. However, the nearly symmetrical broadening around the p pole seen in the 
data is not reproduced by this model. The remaining excess at M > 0.9 GeV/c 2 may 
well be related to the prompt dimuon excess found by NA60 in the intermediate mass 
region 5 . Processes other than 2ir, i.e. An ... could possibly account for the region 
M > 0.9 GeV/c 2 [Ej. 
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3. The anomalous J/ip supression 

In Pb+Pb collisions above a certain centrality threshold an anomalous J/ip suppression 
has been seen by NA50 i- e - suppression mechanisms different from nuclear 

absorption must be invoked to explain the observed 3/ip yield. However, several 
questions raised by this observation have still to be clarified. By studying the J/ip 
production in In+In collisions NA60 investigates the onset of the anomalous J/ip 
behaviour in systems lighter than Pb+Pb. By comparing the suppression pattern 
obtained in different systems as a function of various centrality variables it should 
be possible to single out a scaling variable for the anomalous J/tp suppression. 
Two different analyses have been performed by NA60 to investigate the centrality 
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Figure 4. The 3/ij) over 
DY standard analysis as a 
function of L (see text), 
compared to the results of 
NA50/NA38 IT7I . 
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Figure 5. J/1/1 suppression 
pattern measured in S+U, 
In+In and Pb+Pb, as a 



function of TV, 



part • 



dependence of the J/ip production, corresponding to two different ways of normalizing 
its yield ^Bj- The first technique, the so-called 'standard analysis', is based on a 
normalization to the Drell-Yan events, already widely used in the past by NA38/NA50. 
The study of the ratio between the J/ip and the Drell-Yan cross sections has the 
advantage of being free from systematic errors connected with the efficiency and 
luminosity evaluations. However, the statistical error is large due to the small number 
of high-mass Drell-Yan pairs (a few hundreds); therefore only three centrality bins can 
be defined using the energy measured in the ZDC (-Bzdc)- Figure 0] shows the result 
as a function of the length L of nuclear matter crossed by the charmonium state in 
comparison to the results of NA38/NA50. 

The second technique overcomes the statistical problem by directly studying 
the measured J/ip centrality distribution as a function of Ezdc and correcting the 
inefficiencies due to the reconstruction procedure and the event selection which only 
affect very peripheral events (Ezdc > 15 TeV) and which are very small (< 1%). 
This distribution is compared to the expected pattern if nuclear absorption was the 
only active suppression mechanism. This reference curve is obtained with the Glauber 
model, assuming that 3/ijj production is a hard process, and using cr a b s = 4.18 mb 
for the nuclear absorption cross section of the cc pair ^HJ- In such an analysis, the 



Latest results from NA60 



6 



relative normalization between data and the reference curve is not fixed a priori. We 
have therefore imposed the centrality-averaged ratio between data and reference to be 
equal to the one obtained in the standard analysis. Since the average a j/^/cdy ratio 
has been measured with a 7% statistical error, such an overall error (not shown in the 
figure) applies to the results. 

In Fig. the J /tp suppression pattern measured in In- In collisions is plotted as 
a function of the number of participants -/V par t as determined from Ezdc (see 1201 for 
details). We take into account the experimental resolution of the detector as well as the 
influence of fluctuations in -/V part at fixed impact parameter. The suppression pattern 
is compared to the results obtained for S+U and for Pb+Pb by NA50/NA38 [T7]. The 
observed pattern indicates that a suppression is already present in In+In collisions, 
setting in at ~ 90 participating nucleons. Even if A part seems to be a reasonable 
scaling variable for the anomalous suppression it is fair to say that more accurate 
data for the previously studied collision systems are needed in order to establish a 
firm conclusion on a scaling variable for the suppression mechanism. For comparisons 
as a function of other centrality related variables and for comparisons to theoretical 
predictions see [T%] . 



4. Elliptic flow of charged hadrons and of J/ip 

In non-central heavy-ion collisions collective flow leads to characteristic azimuthal 
correlations between particle momenta and the reaction plane. This is the plane 
defined by the beam direction z and the impact parameter vector b. To quantify the 
anisotropic flow the coefficients of a Fourier expansion of the azimuthal distributions 
with respect to the reaction plane are evaluated. The reaction plane is experimentally 
not directly accessible but it can be determined using the anisotropic flow itself, 
independently for each harmonic n of the anisotropic flow. A detailed overview over 
several methods is given in |2J . The acceptance coverage of NA60 is mainly close to 
midrapidity (y as 3 — 4). Since the directed flow, meaning the first harmonic (n = 1), is 
zero at midrapidity due to symmetry reasons we concentrate on the second harmonic 
(n = 2) only. The method we employ is the so-called event plane method relating 
the azimuthal emission angles <pi of all the charged particles i to the azimuthal angle 
\P of the event plane (an estimate of the reaction plane) by the so-called event flow 
vector Q 

Q x = Qcos(2!?)=^picos(20 i ) , Q y = Q sin(2 V ) = £pj sin(20i) (1) 

i i 

which allows to calculate 

1 f^-iQy\_ 1 f^-i £iMsin(2 



9 = - • tan" 1 ^ = - • tan" 1 rl . (2) 

2 V Qj 2 \ EiKcos(2 •«/»;)/ 

The NA60 experiment has an acceptance which is highly asymmetric in the azimuthal 
angle </>. This strongly affects the distribution of reconstructed event plane angles \P 
as demonstrated by the black circles in Fig. For perfectly symmetric acceptances 
and efficiencies this distributions has to be flat for the average over many events. 
To correct for this asymmetries the event flow vector Q n is recentered in a (pt,y)- 
acceptance matrix by its mean value for each bin of this matrix as averaged over 
many events [22] • This changes Eq. [3] to: 

1 A -i EiPiNn(2-0 i )-(sin(2-0)} pt , y ; 



9 = - • tan" 1 ^ lf ? L ) w /pt,yj (3) 

2 V EiPi[cos(2.0 i )-(cos(2-0)) Pt , y ] / ' 



Latest results from NA60 



7 



x10 3 

2 800 
700 
600 
500 
400 
300 
200 
100 

°-1.5 -1 -0.5 0.5 1 1.5 

Figure 6. Distribution of the reconstructed event plane angles. Black circles: 
without any correction, red squares: after recentering, green triangles: after 
flattening. 

The resulting event plane distribution is shown by the red squares in Fig. [fj] By 
definition the recentering can only remove anisotropics in the same harmonic as the 
one used for the reconstruction of the event plane. To remove remaining anisotropies 
in higher orders a flattening procedure has been applied [2] . This leads to a flat event 
plane distribution as shown by the green triangles in Fig. 
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Figure 7. Azimuthal emis- 
sion pattern of charged par- 
ticles in In+ln at 158 A GeV. 
The line represents a fit to 
the data according to Eg. HI 
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Figure 8. The elliptic 
flow coefficient V2 of charged 
hadrons as a function of pt ■ 



Once the bias in the event plane distribution is removed it is possible to correlate 
the emission angles of single particles with this event plane and to create azimuthal 
emission patterns dN/d(f>, with <f) = '/'lab — 9 being the difference between the emission 
angle </>i a b of the particle as measured in the laboratory and the orientation of the event 
plane >F To avoid autocorrelations the event plane has to be recalculated without 
using the particle track whose emission angle with respect to the event plane is being 
determined. FigurcQshows such an azimuthal emission pattern for charged particles. 
The line denotes the function 

^ ~ 1 + [2 • v' 2 ■ cos(2 • 0)] (4) 
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which was fitted to the data with v' 2 = (cos(2 • <fi)). This coefficient has to be 
corrected for the event plane resolution to obtain the real elliptic flow coefficient 
v 2 = v' 2 /{cos[2 ■ AW}). This resolution has been determined by dividing each 
event in two independend subevents and applying the method described in |23| . 
Depending on the collision ccntrality the event plane resolution varies between 
0.18 < (cos [2 ■AW']) < 0.32. 
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Figure 9. The elliptic 
flow coefficient 112 of charged 
hadrons as a function of pt 
for various collision centrali- 
zes. 
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Figure 10. The elliptic 
flow coefficient V2 of charged 
hadrons as a function of 
the collision ccntrality for 
different ranges in pt ■ 



The resulting elliptic flow coefficients v 2 for all charged particles are displayed as 
a function of pt in Fig. |H1 Please note that the present results are not corrected for 
non-flow correlations due to the HBT effect. Since this mainly affects the results at 
low pt (see e.g. [21]), we only present results for p t > 0.5 GeV/c. Up to now all results 
shown in this section are based on an analysis of approximately 50% of the measured 
statistics. The flow coefficients rise with p t and show a saturation at higher momenta 
as already previously observed in Pb+Pb collisions at the SPS [H] as well as at higher 
energies at RHIC 
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Figure 11. Azimuthal emission patterns of J/t/j in In+ln at 158 A GeV, on 
the left hand side for more central collisions, on the right hand side for more 
peripheral collisions. 
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Using the energy deposited in the ZDC the data can be grouped according to the 
collision centrality. Fig. [5] shows t> 2 for charged particles as a function of p t for four 
centrality bins. The quoted fractions of the total geometrical cross section have been 
calculated within a Glauber model. The trend as a function of pt is similar for all 
collision ccntralities with the values for vi increasing for more peripheral collisions. 
This can be seen in more detail in Fig. 1101 showing vi as a function of the collision 
centrality for three bins in pt . 

Having determined the orientation of the reaction plane with the charged particle 
tracks measured in the vertex tracker we can determine as well the azimuthal angle 
emission pattern of the J/ip with respect to this plane. Figure HP shows the azimuthal 
angle distributions for J/ip, on the left hand side for the more central events (0.5% 

< cr/(jg CO < 20%), on the right hand side for the more peripheral events (20% 

< cr/(jg CO < 83%). In total these distributions contain approximately 12000 J/ip 
which is about 50% of the full statistics available for this type of analysis. The more 
peripheral data seem to indicate a non-isotropic emission pattern but the limited 
statistics docs not allow for a solid conclusion up to now. The analysis of the remaining 
statistics is under way. 
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